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ABSTRACT 

Context. The corona of the Sun is dominated by emission from loop-like structures. When observed in X-ray or extreme ultraviolet 
emission, these million K hot coronal loops show a more or less constant cross section. 

Aims. In this study we show how the interplay of heating, radiative cooling, and heat conduction in an expanding magnetic structure 
can explain the observed constant cross section. 

Methods. We employ a three-dimensional magnetohydrodynamics (3D MHD) model of the corona. The heating of the coronal plasma 
is the result of braiding of the magnetic field lines through footpoint motions and subsequent dissipation of the induced currents. From 
the model we synthesize the coronal emission, which is directly comparable to observations from, e.g., the Atmospheric Imaging 
Assembly on the Solar Dynamics Observatory (AIA/SDO). 

Results. We find that the synthesized observation of a coronal loop seen in the 3D data cube does match actually observed loops 
in count rate and that the cross section is roughly constant, as observed. The magnetic field in the loop is expanding and the plasma 
density is concentrated in this expanding loop; however, the temperature is not constant perpendicular to the plasma loop. The higher 
temperature in the upper outer parts of the loop is so high that this part of the loop is outside the contribution function of the respective 
emission line(s). In effect, the upper part of the plasma loop is not bright and thus the loop actually seen in coronal emission appears 
to have a constant width. 

Conclusions. From this we can conclude that the underlying field-line-braiding heating mechanism provides the proper spatial and 
temporal distribution of the energy input into the corona — at least on the observable scales. 
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1. Introduction 

The corona of the Sun is dominated by loop structures contain- 
ing plasma at temperatures of 10 6 K and above. These coronal 
loops were first revealed by coronagraphic observations in the 
1940s (Bra y et aTT] 1 1 99 1 ). Their highly dynamic nature became 
very clear through extreme ultraviolet (EUV) imaging by the 
Atmospheric Imaging Assembly (AIA; |Lemen et al.||2011|) on 



the recent Solar Dynamics Observatory (SDO; Pesnell et al 
|2011| ). The EUV and X-ray emission from the Sun is closely re- 
lated to the magnetic field in the corona ( [Poletto et al.||1975) and 
most prominently seen above active regions hosting sunspots on 
the solar surface. Stereoscopic EUV observations, together with 
extrapolations of the coronal magnetic field from surface obser- 
vations, have demonstrated that loops roughly follow the mag- 
netic field ( [Feng et aL||2007| ). The formation and appearance of 
these loops in the complex magnetic environment of the corona 
provides a pivotal test for a model of the coronal heating process 
In the corona the energy density of the magnetic field ex- 
ceeds the internal energy density of the plasma, and therefore the 
magnetic field shapes the structures we see in EUV and X-ray 
observations. In contrast to a neutral gas, in an ionized plasma 
the heat conduction is no longer isotropic, but basically parallel 
to the magnetic field. In consequence, neighboring field lines are 
thermally isolated. In addition, the plasma flow has to be parallel 
to the magnetic field, so that one can consider each field line as 
an independent unit. If a part of the corona is now filled with ho t 
plasma, either by chromo spheric evaporation ( Klimc huk] |2006| ), 
injection of heated material from below ( De Pontieu et al.||20lT| ), 
or any other mechanism, only the filled field line will be ob- 
served as a bright coronal loop. In other parts there is still a 



magnetic field, but not enough heating or mass injection to load 
enough plasma onto the field lines. 

Observations show that most of these coronal loops have a 
roughly constant cross section or else expand only little with 
height. This is obvious from the inspection of AIA images 
in the EUV as was hinted at from first subarcesec resolu - 
tion EUV and X-ray images ( [Pere] [1982} |Golub et al.| [T990]), 
and confirmed quantitatively " 



ater (Klimchuk 2000, Watko & 



Klimchuk, 2000). The constant cross section is not simply an ar- 



tifact of image processing, e.g. background subtraction (Lopez 
|Fuentes et al.l|2D08) . 

This constant loop width is in stark contrast to the antici- 
pated magnetic field structure. Extrapolations of the magnetic 
field from the solar surface into the corona show a strong expan- 
sion of the magnetic field. Just as for a magnetic dipole, in a solar 
active region the magnetic field lines should strongly diverge on 
their way into the corona. A quantitative comparison shows that 
the magnetic field lines derived from the extrapolation expand 
by a factor of two or more, while the cospatial EUV loop keep s 
a roughly constant cross section (Lo pez Fuentes et aT) 2006| ). 



If the EUV and X-ray loops follow the field lines as outlined 
above, they should participate in the magnetic expansion, which 
they do not. 

To reconcile the magnetic structure with the EUV and X- 
ray loops numerous suggestions have been made about how 
special magnetic configurations could only show a little expan- 
sion. In force-free, highly sheered magnetic field configurations, 
strong field-aligned currents could prevent the magnetic expan- 
sion ( [McClymont & Mikic 1994). However, the required strong 
currents seem to contradict observations ( [Metcalf et al.| |1994| ). 
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It is also unlikely that additional magnetic twist or a variable 
strong cross-field elongation of the magnetic flux tube defining 
the loop could account for its constant width ( Klim chuk| |2000| ). 
Likewise, for force-free twisted magnetic flux tubes with an in- 
flow of plasma one would not expect a constant cross- section 
loop under solar conditions ( |Petrie||2008| ). Magnetic separators 
connecting magnetic null points can be expected to show less 
expansion with height and thus could host nearly-constant width 
coronal loops (Plowman et al. 2009 ); however, it remains to 
be seen whether coronal loops are always associated with sep- 
arators or with their generalization the (q uasi-) separatrix lay- 
ers within regions of strong magnetic field ( Priest & Demoulin, 
1995). The latter ones are a ssociated with ex panding fans at the 
periphery of active regions ( Schrijver et al. 2010) but have not 
yet been clearly related to loops. 

A different approach is to use an expanding fixed magnetic 
field structure and solve the hydrodynamic problem with a pa- 
rameterized heating function along the field, thereby provid- 
ing the thermal properties to infer the emis sion as observed by 
EUV and X-ray imaging ( |Mok et aL| |2005| ). By optimizing the 
parametrization of the heating function, this procedure can give 
loops of roughly constant width through th e convolution of th e 
resulting temperature and density profiles ( |Mok et aL| |2008| ), 
while it seems that some fine-tuning is required to obtain this 
result. Therefore the logical next step is to account for the full 
interaction of the plasma and the magnetic field and to include a 
more self-consistent treatment of the energy input to sustain the 
hot corona. 



2. 3D MHD coronal model 

In the model we present in this study we incorporate a self- 
consistent time-dependent heating of the plasma in an evolv- 
ing magnetic environment. Horizontal granular motions in the 
photosphere braid the magnetic field lines which leads to cur- 
rents in the upper atmosphere that subsequently heat the plasma 
through Ohmic dissipation. This concept follows the nanoflare 
model ( [Parker |1972[ |1983|) and is similar to the flux-tube tec 



lated. To solve the problem num erically we use the Pencil codqj 
( [Brandenburg & Dobler[[2002| ). 

The computational domain covers 50x50 Mm 2 in the hori- 
zontal and 30 Mm in the vertical direction. The grid spacing is 
230 km in the vertical and 390 km in the horizontal direction. For 



the initial condition we use the same magnetogram as Bing ert & 
Peter ( 2011) and fill the volume with a magnetic field based on 
a potential field extrapolation. For the initial temperature we as- 
sume a simple solar-like average stratification and calculate the 
initial density from hydrostatic equilibrium. We then force the 
lower boundary by applying a granulation-type horizontal mo- 
tion as described by Bingert & Peter (2011). This has the same 
statistical properties as real granulation. At the bottom the mag- 
netic field is frozen to the plasma motions, which leads to the 
braiding of the field lines. In the horizontal direction we employ 
periodic boundary conditions. 

As in other models of the same type (e.g. IGudiksen & 



Nordlund , 2002] |Hansteen et aT) [20T0l |Bingert & Peter ||20 11 



Ohmic heating dissipates the induced currents j and heats the 
corona. This is through a term r] /io j 2 in the energy equation 
with the magnetic permeability /io and a constant magnetic re- 
sistivity ?7=10 10 m 2 /s. As in the aforementioned models we use 
a value for 77 so that the magnetic Reynolds number when us- 
ing the grid spacing for the length scale is on the order of unity, 
which requires rj to be much larger than what follows from clas- 
sical transport theory (e.g. |BoycT& Sanderson, 2003). This en- 
sures that current sheets have a thickness of at least the grid spac- 
ing, so that the numerical scheme can resolve them. On the real 
Sun, for much lower values of 77, much thinner current sheets 
would form. Still, as long as the Poynting flux into the corona is 
the same, the dissipated ener gy should remain the same, because 
the currents are then higher ( Hendrix et al. 1996[ |Galsgaard &| 
Nordlund] 1 1996| ). Of course when going to smaller and smaller 
length scales, other physical processes come into play, especially 
when going below the electron mean free path (several 10 km) 
or the electron gyro radius (several m). In this study we are in- 
terested in the observable consequences on currently resolvable 
scales, so one should consider the Ohmic heating term 77 /io j 2 as 
a parameterization of the true heating mechanism. 



tiestet ai.||2002] ). Numerical experiments using 3. Synthesized and real EUV images 



tonics model 

three-dimensional magnetohydrodynamics (MHD) models show 



that this process leads to a loop-dominated corona (Gudiksen 
& Nordlund, 2002, 2005a b) and that many derived observable 



quantities match observed average properties ( 


Peter etal.||2004| 


2006J Hansteen et al.| 2010( Zacharias et al. 


201 1}, including 



the transition region Doppler shifts or the differential emission 
measure. 

The details of our numerical model have been described 
before (Binger t~& Peter| |2011| ). Basically, we solve the time- 
dependent 3D MHD problem, i.e. the induction equation along 
with the mass conservation, momentum, and energy balance 
from the cool solar surface into the hot corona above an ac- 
tive region. The energy equation includes Spitzer heat conduc- 
tion along the magnetic field, optically thin radiative losses and 
Ohmic heating. Only this ensures a proper description of the 
coronal pressure, which is prerequisite to synthesizing EUV 
emission as would be expected for, say, AIA observations. For 
numerical stability also a (small) heat transfer perpendicular to 
the magnetic field is needed. In the coronal part of the compu- 
tational domain the perpendicular heat conduction is about 100 
or more times weaker than the parallel conduction, which en- 
sures that neighboring magnetic structures are still thermally iso- 



In Fig.[l]we display the synthesized emission in the AIA 171 A 
channel from the computational domain showing plasma at just 
below 10 6 K. This image was synthesized using the AIA tem- 
perature response functions ( [Boerner et al.| |2011[ |Peter et aL 
|20T2] ). 

A coronal loop sticks out, connecting the periphery of the 
two main polarities (sunspots) in the photosphere. This shows 
only a snapshot from the numerical experiment, and the struc- 
tures quickly evolve in time — some 30 minutes later other 
structures close-by will dominate. The image also changes with 
the viewing angle, of course, as a result of the line-of-sight inte- 
gration of the optically thin emission. The figures shown in this 
paper show snapshots from the simulation, all about one hour 
solar time from the start of the simulation. The temporal evolu- 
tion for some 50 minutes after this snapshot is available in the 
online-edition. 

For comparison in Fig. [2] we show an actual observation of an 
active region using AIA. A system of loops connecting the main 
polarities of the active region is dominating the field-of-view. 
The numerical experiment presented in this manuscript was not 



http://code.google.eom/p/pencil-code/ 
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Fig. 1. View of the modeled corona above an active region as 
would be seen by AIA/SDO in the 171 A channel dominated by 
emission from below about 10 6 K. The rectangle indicates the 
boundaries of the computational domain (50x50x30 Mm) at 
the coronal base. The most prominent loop connects the periph- 
ery of the two magnetic concentrations at the surface (not shown 
here) and has a roughly semicircular shape with about 30 Mm 
length. In the vicinity of the active region, hazy emission can 
be seen that is associated with the diffuse background corona in 
the quiet Sun. The white dashed line indicates the position of a 
loop investigated further in Fig. [3] The temporal evolution over 
50 min is shown in a movie available in the online edition and 
at http://www.mps.mpg.de/data/outgoing/peter/papers/2012- 



loops/aa 1 947 3 -peter-fig 1 .mp4 



set up to match the specific observation shown here. A movie 
showing the temporal evolution is available in the online edition. 

The count rates of the coronal images synthesized from the 
3D MHD active region model roughly match the typical count 
rates actually observed in active regions with AIA. The count 
rates in the synthesized loop are typically a factor of two lower 
than in the loops of the real active region (see color scales in 
Figs.[T]and|2]and the profile in Fig.[3]). This difference in count 
rate relates to a difference in density of about a factor of 1.4. 
This is small considering the scatter of active region brightness 
and the fact that we did not specifically designed the model to re- 
produce this particular active region. This shows that the model 
loop has roughly the same amount of emitting material as the 
real Sun, therefore the thermal structure (density and tempera- 
ture) found in our model should be similar to the Sun. Besides 
the well defined loop we also see other structures, some with 
higher contrast, some appearing quite diffuse, similar to what is 
found in real observations. Of course, the actual active region 
shows a finer structure, which reflects the limitations of the cur- 
rent MHD models. 

The temporal evolution of the loop in the synthesized corona 
and the real observed loops (see online material with Figs.[T]and 
[2]) is similar in the sense that in both cases the loops have a lim- 
ited lifetime of some 20 to 30 minutes. However, the real ob- 
servations show a higher variability, especially on smaller scales 
that are not described by our model. This again reflects the lim- 
itations of the model. But still, the model and the observations 
have the common property of a lifetime of a larger structure of 
about one half hour, which is consistent with the coronal (radia- 
tive) cooling time. Further detailed work is needed, especially 
for more complex and longer numerical experiments, to draw 



Fig. 2. Real observation of an active region. This snap- 
shot shows an active region observed with AIA in the 
171 A bandpass on Aug 08 2010 at 01:13:49 UTC. The 
field-of-view covers 150" x 90" centered at solar (X,Y) 
« (185", 310"). South is top. The white dashed line indi- 
cates the position of a loop investigated further in Fig. [3] 
The temporal evolution over 75 min (from 00:36 to 01:51) 
is shown in a movie available in the on-line edition and 
at http://www.mps.mpg.de/data/outgoing/peter/papers/20 1 2- 
loops/aa 1 947 3 -peter-fig2 .mp4 
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Fig. 3. Intensity variation along selected coronal loops. The bot- 
tom curve shows the variation along the loop synthesized from 
the 3D MHD model outlined in Fig. [I] The top curve shows the 
same for the observed loop outlined in Fig. [2] The ordinate is 
plotted logarithmically. 



further conclusions on how well the models reproduce the vari- 
ability seen in coronal observations. 

The intensity distribution along the loop synthesized from 
the model (and highlighted in Fig.[T]) is shown in Fig. [3] The foot- 
points are slightly brighter compared to the loop apex. This can 
be compared to the intensity variation for one particular loop in- 
dicated in Fig.|2] which is overplotted in Fig. [3] The ratio of the 
count rate for the observed and synthesized loop is about a fac- 
tor of two. Thus the variation along the synthesized loop roughly 
matches observations; however, as the loops in the model and 
the observati on change in time, this can only be a snapshot. Mok 
|et al.| (|2008|) report in their loop study that "an important defi- 



ciency of [their] model, as well as other models (Schrijver et al.| 


2004, 


Brooks & Warren 


2006 ; Warren & Winebarger 


2006), is 
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(a) synthesized AIA 171 A 




the excessive EUV emission from the footpoints." In the model 
presented in this study we have not found these overly bright 
footpoints, but the ratio of loop apex to footpoint is roughly 
comparable to observations. This is not only true for the snap- 
shot shown in Figs.p]and[3| but it also holds for other times (see 
animation attached to Fig]j]m the onlin e edition). One re ason 
for this difference between the model of |Mok et al.| ( [2008] ) and 
our model could be that theirs has a stronger concentration of 
the heating rate (per particle) towards the loop footpoints, which 
could lead to an increased radiative output near the loop feet (see 
discussion in Appendix[A|). 

4. Loop cross section 

Already in Fig. [T] the loop looks like is has a more or less con- 
stant cross section. This becomes more evident when integrating 
through the box horizontally (along a line-of- sight perpendicular 
to the loop plane). This is shown in Fig.|4]where we also applied 
the point-spread-function (PSF) for the AIA 171 A channejjfor 
a better comparison with real AIA images. The image is dis- 
played with the actual AIA pixel size and after background sub- 
traction (as also used in observations to enhance the contrast). 
From this image one would roughly estimate the loop width by 
eye to be some four to five AIA pixels wide, corresponding to 
some two to three Mm. For a roughly 30 Mm long loop (such as 
the one in Fig. [4]), this is consistent with observations and with 
an order-of-magnitude estimation for the loop wid th based on 
the field-line-braiding mechanism ( Schrijver, 2007 ). 

To quantify the width of the coronal loop and its expansion 
we investigate the full width at half maximum (fwhm) of the 
emission as a function of arc length along the loop. In Fig. [5^ 
we show the same as in Fig.[4](now without the PSF applied, i.e. 
with the full model resolution), together with the magnetic field 
line along the spine of the loop. For five selected cross sections 
we show the (normalized) intensity profiles in panel b. The (lin- 
ear) expansion of the loop is defined as the FWHM normalized 
to the fwhm at the loop top and is marked for the selected pro- 
files in panel c by the crosses. The red solid line connecting the 
crosses shows the expansion all along the loop. For comparison 
the dashed line shows the expansion if one applies the AIA PSF. 
Either way, the linear expansion is small, of the order of 20% 
and less. In contrast, the magnetic expansion is much stronger. 
Because the magnetic flux is conserved, the area expansion is 
proportional to the inverse of the magnetic field strength \B\. 
Thus the linear (or radial) magnetic expansion is proportional 
to 1/ which is shown by the black solid line in Fig. [5J:. 
Clearly, the (linear) magnetic expansion is well above a factor of 
two from the coronal base to the apex. 

From this we can conclude that we see a clear expansion of 
the magnetic field, similar to a potential field extrapolation. This 
is to be expected for the flux braiding mechanism: the magnetic 
resistivity will prevent the build-up of too large twist through 
the footpoint motions and thus the magnetic field stays close to 
poten t ial ([Ga lsgaard & Nordlund, 1996; Nordlund & Galsgaard 



|1997[ |Schrijver| |2007| ). Nonetheless, the resulting synthesizec 
EUV loop shows a constant cross section. 

5. Discussion 

To understand the nature of the constant cross section of the 
loop we investigate the plasma properties in a 2D cut through 
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The point spread function (PSF) for AIA is available in SolarSoft 
( http ://w ww.lmsal.com/solarsoft/ ) . 



Fig. 5. Analysis of the width of the synthesized coronal loop. 
Panel (a) displays the emission (same as Fig.[4] now without the 
point spread function, PSF). Overplotted are the magnetic field 
line along the loop spine (red dashed) and selected cross sec- 
tions. Panel (b) shows the (normalized) cross-sectional intensity 
profiles along the lines in panel (a) with matching color coding. 
Panel (c) shows the width of the loop normalized to the width at 
the apex, which defines the (linear) expansion of the loop. This 
expansion with (dashed) and without (solid) the PSF applied to 
the synthesied image is plotted in red. Overplotted is the mag- 
netic expansion (black; cx|I?| -1 / 2 ) along the spine field line. The 
crosses indicate the expansion at the cross-sectional profiles in 
panel (a). 

the loop plane which is shown in Fig. [6] While panels a, b, and d 
show the properties at a given time in the simulation, the heating 
rate in panel c shows the time average over ten minutes before 
the snapshots in the other panels. The heating rate per particle 
is concentrated towards the footpoints, consistent with earlier 
investigations. The increased heating rate, roughly confined to 
the magnetic loop defined by the field lines, leads to enhanced 
heat conduction back to the Sun, resulting in evaporation and fill- 
ing the heated structure with mass, as visible in Fig. [6^. The en- 
hanced density defines a plasma loop, which is expanding along 
the magnetic field. 

The temperature structure, however, is not aligned with the 
magnetic field (Fig. [6})). Each field line is basically thermally 
isolated from the neighboring field lines and along each field line 
the temperature rises towards the apex (similar to what is found 
in ID models). Because the heating rate is not the same for each 
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Fig. 4. View of the synthesized coronal loop in 
AIA 171 A seen in Fig.[I]when horizontally in- 
tegrating through the box along a line-of- sight 
perpendicular to the loop plane. The AIA 171 A 
point spread function (PSF) was applied and 
shows the image with the same pixel size as 
AIA (0.6 arcsec corresponding to 435 km on 
the Sun). As in observations, a background sub- 
traction was applied to enhance the contrast. 
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Fig. 6. Two-dimensional cut through the loop plane. Panels (a) and (b) show the density and temperaure at a snapshot at the same 
time as Figs.[T]to|5] Overplotted are the magnetic field lines (dashed) defining the magnetic loop, which roughly coincides with the 
plasma loop defined by the density structure. Panel (c) displays the heating rate per particle averaged over 10 minutes before the 
time of the other snapshots. This also roughly coincides with the magnetic loop. Panel (d) shows the synthesized emission in the 
AIA 171 A channel defining the eLoop outlining the EUV emission. This is smaller than the magnetic loop and the plasma loop and 
has roughly constant width. The lines labeled [a] and [b] indicate the cuts shown in Fig. [7] 
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Fig. 7. Profiles of AIA 171 A temperature response, temperature, 
density, and AIA 171 A emission along the cuts [a] and [b] in 
Fig. [6] The horizontal bars show the full width at 80% of the 
maximum value (numbers giving the values in Mm). 



field line also the temperature on different field lines will not be 
the same, just as if a collection of field lines would be treated 



in individual ID models. As in scaling laws (e.g. Rosneretal 
|1978| ) longer field lines will reach higher temperatures (for the 
same heating rate). This leads to a temperature gradient across 
the loop structure as visible in Fig.[6]3 and[7]with the temperature 
increasing "outwards". 

The AIA response function for the 171 A channel has its 
maximum around 0.9 MK, i.e. in the reddish region in Fig.JSJ). In 
the upper part of the magnetic loop the temperature is too hot for 
the AIA 171 A channel and despite the significant density, there 
the AIA emission is only small. Consequently the AIA emission 
is confined to the lower part of the magnetic loop (Fig.[6]i). The 
resulting loop defined by the emission, the "eLoop", shows a 
more or less constant cross section. 

To quantify the formation of the eLoop, we show in Fig.|7] 
cuts across the loop from Fig. [6] In the vertical cut [a] the tem- 
perature profile (across the field) restricts the region, which con- 
tributes to the AIA emission to the lowermost part of the density 



peak; consequently, the emission profile is much narrower than 
the density profile. For the inclined cut [b], the AIA contribution 
is quite flat across the density peak and thus the emission profile 
basically follows the density squared. Between the two cuts [a] 
and [b], the magnetic field strength on the loop spine drops from 
about 30 G to 16.5 G, which corresponds to a (linear) expansion 
of the magnetic field by some 35%. This is basically the same 
for the expansion of the density structure, i.e. the plasma loop 
(from 3.6 Mm to 4.9 Mm). However, the expansion of the emis- 
sion defining the eLoop is only some 3% from about 2.5 Mm to 
2.6 Mm. The width of the plasma loop and the eLoop has been 
determined by the full width at 80%. In contrast to the analysis 
in Fig. [5] no background subtraction was applied here leading to 
less contrast between the region in and next to the loop. 

So far we have presented results only for the AIA 171 A 
channel, because this one shows the least contamination by cool 



material of the coronal AIA channels (Martmez-Sykora et al. 



|2011| ). If we repeat the analysis for the AIA channels that show 
cooler and hotter plasma, e.g. the 13lA or 193 A channels, 
we also see a loop at constant width (see Fig. [8}. Their cross- 
sectional profiles are more or less Gaussian with a width of 
some 2" to 3", c onsistent with observations of, e.g., Wat ko & 
Klimchuk (2000). They have a small spatial offset to the 171 A 



loop by some 1/4" to 1/2", but it is still overlapping. The small 
offset in location of the loops in the passbands is consistent with 
observations that often show almost cospatial emission of loo ps 
in the different AIA channels ( Aschwanden & Boe rner[|2011| ). 

The fundamental reason for the constant cross section found 
in this study is the presence of a (small) temperature gradient 
across the structure of enhanced density. In combination with 
the narrow contribution function to the EUV lines (and thus the 
EUV passbands of AIA), this cuts off part of the structure with 
enhanced density, as illustrated in Figs. [6] and [7] Therefore our 
explanation applies not only to the coronal bandpasses of AIA 
discussed here, but also to the previous TRACE and EIT/SOHO 
data. It even applies to single line rasters with EUV spectrom- 
eters such as EIS/Hinode, because the (temperature) contribu- 
tion functions of individual coronal lines typically have widths 
of some 0.3 to 0.4 in logT[K], which is comparable to the 
widths of the peaks of the AIA temperature response functions 
(cf. |Boerner et al.||2011| ). In principle this argumentation would 
also apply to the formation of transition region loops (seen e.g. 
in C 4), but the model presented here does not provide the reso- 
lution needed to properly resolve these cool structures. 

At this point we cannot say much about the mechanism that 
leads to the constant cross- section of loops seen in X-rays. A 
comparison of synthesized X-ray images to observations would 
be of particular interest, because the X-ray filters cover a much 
wider range of temperatures than the EUV filters. However, 
the contribution functions of the filters in current and past X- 
ray instrumen ts mostly peak well above in logT [K]=6.5 (e.g. 
Takeda] |2011| ). This is beyond the peak temperature we reach 
in the model presented here. Therefore future models reaching 
higher temperatures will have to show whether the explanation 
provided here for the EUV images will also hold for the X-ray 
images. 



6. Conclusions 

We have presented a self-consistent 3D magnetohydrodynamics 
model of the solar corona, which for the first time gives a de- 
tailed explanation of loops in the corona that ubiquitously show 
a nearly constant cross section in EUV observations. Through 
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Fig. 8. Appearance of loop in different passbands of AIA. Panels 
(a) to (c) display the loop seen from the side in AIA 131, 
171, and 211 A. The peak contribution of these bands is at 0.6, 
0.8, and 1.9 MK, respectively. The images are shown with the 
AIA pixel size, and with the AIA PSF applied and background- 
subtracted. The horizontal dashed line is plotted for reference 
only. The bottom panel (d) shows the variation across the loop 
near its apex as indicated by the vertical line in the top pan- 
els (ordinate plotted logarithmically). The vertical lines in panel 
(d) indicate the position of maximum emission in the respective 
band. 



this we solve a problem that has remained enigmatic since the 
first coronal loop observations. 

We find that the cross- sectional intensity profiles are close to 
Gaussians with widths of a few a rcsec, just as found in observa- 
tions ( |Watko & Klimchuk, 2000). While the loop seen in emis- 
sion has h ardly any expansion, again as in observations (Watko 
& Klimchuk, 2000), the plasma fills the magnetic structure; i.e., 



the plasma loop closely follows the expanding magnetic field. 
The constant width of the loop is then due to the variable con- 
tribution of the plasma at different temperatures. In a manner 
of speaking, the contribution function (depending on tempera- 
ture) cuts off parts of the expanding plasma loop, and the re- 
maining loop seen in emission then has a more or less con- 
stant cross section. From these results it is also clear that a ID 
model for a coronal loop would not be able to capture the mecha- 
nism of the formation of a constant-cross- section loop as outline 
here. Furthermore, the loop width obtained through our model 
is consistent with EUV loop observations and a simple order-of- 



magnitude model for the field-line-braiding process (Schrijver 
|2057> . 

All together, this provides evidence that the process of field- 
line braiding with subsequent dissipation of the induced cur- 
rents, which is at the heart of our model, provides the proper 
spatial and temporal distribution of energy input to heat and 
drive the dynamic corona on observable scales. Of course, this 
does not exclude other mechanisms, but shows that the 3D mod- 
els based on field-line braiding passed yet another observational 
test. 
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Fig. 9. Heating rates per particle. The 2D histogram displays the 
distribution of the Ohmic heating rate per particle, with the hor- 
izontal averages shown as a solid line. Overplotted as a dotted- 
dashed line is the horizontally averaged parameterization of the 
heating rate as used by Mo k et al.| p008). 



Appendix A: Heating rate in the 3D MHD model 

At the end of Sect.[3] we mentioned that one possible cause for 
the overly bright footpoint regions of the loops found in previous 
models (Mo k et aT| |2008| and references therein) might be the 
strong er concentration o f the heating rate towards the footpoints 
in the |Mok et ah] ( |2008| ) model. To substantiate this, we give a 
short compa rison of the heating rates in our model and in Mok 
[etaLl (2055) . 

In our model the energy input is due to Ohmic heating 
through the dissipation of currents j induced by field line braid- 
ing with a constant magnetic resistivity r] at a rate 



In Fig. [9] we show the 2D histogram of the Ohmic heating rate 
per particle (volumetric heating rate divided by the particle den- 
sity) in our model as a function of height in the atmosphere. For 
this we evaluated i^ohm at each grid point for one snapshot. The 
picture does not change (significantly) over time. Overplotted 
as a thick solid line is the horizontally averaged Ohmic heating 
per particle. As emphasized in previous studies, the heating per 
particle shows a pe ak at the base of th e corona. 

In their model Mok et al. (2008) use a parameterization of 
the heating rate through the local magnetic field strength B and 
density p, 

^Mok oc B 2 5 



-0.5 



Again, we evaluate i^Mok at each grid point and overplot the hor- 
izontal averages of i^Mok per particle as a function of height in 
Fig-El 

It is clear from F ig. [9] that the heating rate per particle as used 
by |Mok et al.| ( |2008) shows a stronger concentration towards the 
base of the corona. This could be part of an explanation of the 
too strong e mission near the loop footpoints as found by Mok 
|eTaLl ( [20n8l ). 
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